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We discuss the non-zero baryon number formulation of QCD in the quenched limit at finite temperature. This 
describes the thermodynamics of gluons in the background of static quark sources. Although a sign problem 
remains in this theory, our simulation results show that it can be handled quite well numerically. The transition 
region gets shifted to smaller temperatures and the transition region broadens with increasing baryon number. 
Although the action is in our formulation explicitly Z{3) symmetric the Polyakov loop expectation value becomes 
non-zero already in the low temperature phase and the heavy quark potential gets screened at non-vanishing 
number density already this phase. 



1. Introduction 

An important aim of lattice QCD is the under- 
standing of the QCD phase diagram and its de- 
pendence on the temperature T and the baryon 
density ns- Especially the region of non-zero den- 
sity is important, as it describes the behaviour 
of dense matter created in heavy ion collisions 
and plays an important role in the cosmological 
context. Due to the well known problem of the 
complex fermion determinant when a non-zero 
chemical potential /i is introduced |^,^ , only qual- 
itative features of the phase digram at non-zero 
density can be understood in terms of models and 
approximations. 

Introducing a chemical potential /x [|| leads to the 
grand canonical partition function of finite den- 
sity QCD. An alternative formulation is given in 
terms of the canonical partition function at fixed 
non-zero baryon number Q . This is achieved by 
introducing an imaginary chemical potential [^,^ 
in the grand canonical partition function and per- 
forming a Fourier integration to project onto the 
canonical partition function for a given sector of 
fixed baryon number W 



Z{B,T,V) 



1 

2^ 



'^'l'Z{i<j>,T,V). (1) 



This formulation still suffers from a sign problem, 
but leads to a quite natural and useful formula- 
tion of the quenched limit of QCD at non-zero 
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baryon number density, where the sign problem 
can be handled quite well 

2. The quenched limit of QCD at non-zero 
density 

The static limit of QCD at non-zero chemical 
potential /i in the grand canonical approach has 
been formulated in |^ and j^. The numerical 
results indicate that in this case the first order 
deconfinement transition of the SU{3) gauge the- 
ory turns into a crossover for arbitrarily small, 
non-zero values of the chemical potential j|] . The 
quenched partition function in the canonical ap- 
proach can be written as 

Z{B,T,V)= fUdU^^Jse-^'' (2) 

where the constraint on the baryon number is en- 
coded in the function fs which is a function of 
Polyakov loops, B counts the number of quarks, 
i.e. B/5 is the baryon number. For B — 3 and 
one flavour of Wilson fermions fs is, for instance, 
given by 



+F(12 + ^[L3,o]-2[Li,i])) 



(3) 



with [L,,,] = T/-i^-(Trij)^(Tri|p. 
For a more detailed description of fs see 



1 

0.8 
0.6 
0.4 
0.2 




m 



B=6 8'^x2 * 
B=12 ^fx2 



B=12 8 x2 



4.854.94.95 5 5.055.1 5.155.25.255.3 



0.5 
0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




<|L|>v 



IT" 



■ B=0 8''x2 + 
■ i * B=0 10-^x2 • 
..■ B B=6 8h2 ^* 

B=12 10^x2^ 
'aiB* * B=12 8''x2 ■ 

- + + + + + + + + ^ 



4.854.94.95 5 5.055.15.155.25.255.3 









t| 16^x4 B=6» 




I 12^x4 B=6" 








P 



5.55 5.6 5.65 5.7 5.75 5.8 



Figure 1. (sgn(Rc/B)) for _B = 6 and 12 calcu- 
lated on lattices of size x N^. 

Sa is the gluonic action, which is Z{S) symmet- 
ric. The partition function Z{B,T,V), is non- 
zero only if _B is a multiple of 3, because fs is 
invariant under Z{3) transformations only if B 
is a multiple of 3. In general it changes by a 
factor e2iriB/3 ^jjjjgp ^ global Z{3) transforma- 
tion of time-like link variables. fs is still a 
complex function, but upon integration over the 
gauge fields the imaginary part of the partition 
function vanishes. The remaining sign problem 
can be handled by using the absolute value of 
Re/s and including the sign in the calculation of 
observables (I). 

Our simulations are performed on x Nr lat- 
tices with = 8,10,12,16 and Nr = 2,4 us- 
ing the standard Wilson action and one flavour 
of Wilson fermions with quark number values of 
5 = 6 and 12 at fixed ub/T^ = {l/3)B{Nr/N^f. 
Fig. 1 shows the average sign (sgn(Re/B)) as a 
function of the coupling /3. For large values of 
the temperature the sign is almost always posi- 
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Figure 2. Polyakov loop expectation value 

for different values of B and lattices of size x 

Nr. 

five, but also for the smallest temperature in our 
analysis the sign can be well determined. It de- 
pends on the spatial volume N^ but varies little 
with B. The Polyakov loop expectation values in 
Fig. 2 show a clear signal for a first order transi- 
tion for the B ~ case, while for all _B > the 
transition is continuous. The transition region 
is shifted towards smaller /3- values and it broad- 
ens with increasing B. Note that by changing 
the gauge coupling /3 we vary the lattice cut-off 
and through this also the baryon number density 
continuously. The broadening of the transition 
region may indicate the presence of a region of 
coexisting phases. The Polyakov loop susceptibil- 
ity (Fig. 3) reflects the existence of a transition 
region that becomes broader with increasing ub, 
but does not show indications for a discontinuity. 
The Polyakov loop expectation value becomes 
non-zero already in the low temperature phase. 
This indicates that the heavy quark potential 
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Figure 3. Polyakov loop susceptibility xl for dif- 
ferent values of B and lattices of size x Nr- 

stays finite at large distances. We validate this 
by calculating the potential using Polyakov loop 
correlations (Fig. 4). For zero baryon number it 
shows the usual linearly rising behaviour for the 
quenched case. For B = 6 the potential stays fi- 
nite at large distances due to the screening of the 
static quark anti-quark sources by already present 
static quarks. This behaviour is comparable to 
heavy quark potentials in full QCD pO| |. 

3. Conclusions 

We have analyzed the quenched limit of QCD 
at non-zero baryon number. The sign problem in 
this theory can be handled quite well numerically. 
We find indications for a region of coexisting 
phases, which broadens with increasing baryon 
number density and is shifted towards smaller 
temperatures. Further analyses are needed to 
see if this is a signal of the existence of a first 
order phase transition or a smooth crossover at 
non-zero density. We also see evidence that the 
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Figure 4. Heavy quark potential for f3 = 5.62 and 
B = and 6. 

heavy quark potential for non-zero baryon den- 
sity stays finite for large distances already in the 
hadronic phase. The potential gets screened by 
the static quarks that induce the non-vanishing 
density. This will have a direct influence on heavy 
quark bound states at high density. 
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